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In this short review, we discuss recent advances in exact solutions of models based on a one- 
dimensional (ID) Coulomb gas by means of field-theoretic functional integral methods. The exact 
solutions can be used to assess the accuracy of various approximations such as the weak coupling 
Poisson-Boltzmann theory as well as the strong coupling theory of Coulomb gases. We consider three 
different ID models: the Coulomb fluid configuration in the case of the soap film model consisting 
of positively and negatively charged particles between adsorbing boundaries, counter ions between 
two charged surfaces, and an ionic liquid lattice capacitor with positively and negatively charged 
particles on a lattice between one positive and one negative bounding surface. 

I. INTRODUCTION 

Field-theoretic functional integral methods can be used to study exact solutions of models based on a one- 
dimensional (ID) Coulomb gas with charged boundaries. In ID, exactly solvable Coulomb gas models can be 
then used as a testbed for assessing the accuracy of various approximations: the weak coupling expansion, Poisson- 
Boltzmann/mean- field equations, and the strong coupling expansion pQ. We review these approximations in the 
context of three ID Coulomb gas systems and remark on whether or not they fail to predict important effects present 
in the exact solution. 

Some physical properties of the ID system can be applicable at least qualitatively for dimensions d > 1 and can 
help us to understand whether pertaining approximation methods are reliable or not. In particular, our analysis gives 
insight into systems such as an array of charged smectic layers or lipid multilayers, and ionic liquids near charged 
interfaces, treated as effectively ID systems. An important aspect of these endeavours is that we can test and develop 
the analysis and especially numerical methods that can then be tentatively applied also for d > 1. 
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FIG. 1: Top: The ID Coulomb fluid configuration in the case of the soap film model (positively and negatively charged 
particles and adsorbing surfaces), middle: counterions between charged surfaces (positively charged particles and oppositely 
charged bounding surfaces) and bottom: an ionic liquid lattice capacitor (positively and negatively charged particles on a 
lattice with one positive and one negative bounding surface). Positively (negatively) charged particles are shown schematically 
as light (dark) gray spheres. 



II. THEORETICAL METHODS 

The method of functional integrals applied to Coulomb gas systems has been developed over many years [2HS]. 
In any dimension this approach allows for both strong and weak coupling to be studied explicitly, but specifically 
in ID the functional integral representation can be applied using a variety of methods to obtain exact solutions to 
a number of models which are generally characterized by a Coulomb gas of ions of possibly non-zero size confined 
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between boundaries with properties that allow their potential or charge to be determined either dynamically or as an 
external field condition. Three varieties of a ID Coulomb gas model discussed below are presented in Fig. [I] 

The functional integral representation of the Coulomb gas partition function allows us to formulate two effective 
solution techniques. The Schrddinger kernel technique is applicable in all dimensions and has been used to analyze 
a number of models [6]. In ID it corresponds to solving the Schrodinger equation [2 which is in principle exact. 
In d > 1 the Schrodinger kernel field theoretic representation of the partition function is derived, often using a 
Hubbard- Str at onovich transformation, and is analyzed by perturbative and graphical methods. For d > 1 this 
approach does require that a preferred co-ordinate can be designated as the Euclidean time and so the approach is 
limited to symmetrically layered systems [6 . The transfer matrix and Fourier methods technique is an alternative to 
the Schrodinger kernel approach. Though it is more general, it is only practical in ID. Its implementation exploits 
periodicity in the (imaginary) electrostatic potential <j) which also restricts its general applicability. Examples of this 
technique in ID are the counterion gas and lattice ionic liquids. 

The actual formulation of the functional integral method relies on the action for the full QED of a general system 
that is then reduced to the electrostatic action proper. The relevant electrostatic Lagrangian is then 

CM = \ej dx(V^(x)) 2 - e^g^(xi) - J dxp e (x)^(x). 

i 

Here qi is the charge of the z-th ion at position x^ and p e ( x ) is the external charge distribution. The partition function is 
obtained by tracing the Boltzmann weight of the above Lagrangian over the electrostatic field [ip] . Tracing furthermore 
over ion positions, changing the axis of functional integration if; = i<p and introducing fugacity p- = p+ = p by the 
Gibbs technique, the partition function for monovalent ions (with qi = ±1) assumes the form 

Z = J d[4>]exp(S (</>)) , (1) 

with the "field action" 

S((j)) = -^e J dx(V0(x)) 2 + 2/i J dx cos(e/?0(x)) , (2) 
where f3 = 1/(&bT). The charge density operator is then given by 



P = p-^\ogZ(p) p = 2/i(cos(e/30)), (3) 



where (• • • ) stands for the <fi average. 



III. BILAYER SOAP FILM IN IONIC SOLUTION 



Because hydrophobic heads of the surfactant molecules preferentially migrate to the surfaces charging them up 
dynamically, the configuration of the bilayer soap film consists of two planar (surfactant) surfaces separated by a 
distance L confining a solution of a symmetric electrolyte. 

We calculate the surface charge, the density profile of electrolyte near the interfaces, and the disjoining pressure P 
as a function of the thickness L of the soap film, defined as 

P-Pmm-Pbuik--^^ ^ 

i.e. the difference between the film and bulk pressures. Here J is the grand-canonical partition function/unit area. 
An important phenomenon to predict is the first-order collapse transition of the film to a Newton black film expected 
as the electrostatic coupling in the film is increased. 

We model this system by a Coulomb gas confined to z G [0,L], schematically presented in Fig. [l] (top), with 
potentials on the boundaries that acc ount for the hydrophillic nature of the head group of the surfactant molecule. 
The Debye length is given by Id = \/ eksT /2pe 2 , and the Bjerrum length in ID by Ib = 2/c#T£/e 2 . Perturbation 
theory is an expansion in the coupling parameter g = Id/Ib- 

We use the partition function described earlier but now includes surface free energy f((j)) to model the surface 
potentials, which are attractive for the negatively charged hydrophillic surfactant head groups whose surface density 
is denoted by p-((j)): 
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where A controls the potential strength. To simplify the notation we scale the variables: <fi — »• e/30, x — x/#. The 
charge density operators for ± charges are then given by the Boltzmann weights p±(<fi) = er^^ . The ID partition 
function then becomes 

where 

K{<j)Q,<l) x \x) = J V(j)(x) exp y dx f C{(j)(x f )) 
is the Schrodinger kernel for evolution in the "Euclidean time" x: 

*(0,a;) = j dcj>'K{4>, ( j> l ;x)f{4>'). 
It satisfies the Schrodinger (Feynman-Kac) equation 

Pe*dx ' + 2^ C0S ^' 

with Z(#) = l/(cos(0)). The above equation is also known as the Mathieu equation, and the harmonic term gives the 
Debye length in units of Ib- Z(g) = 2/i/p is the renormalization that relates the fugacity to the observable charge 
density and is given by Eq. We now consider the solution in various limiting regimes. 



A. Large L: bulk pressure 

Strong coupling (SC) g —> oo: The Mathieu ground state dominates in this regime and so we can use the 
Schrodinger perturbation theory for the ground-state energy of H. The result, derived originally in [2], is 
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The leading term is the free gas term but for density p/2, which therefore signals the onset of the dimerization process, 
i.e., the Bjerrum pair formation of positive and negative mobile charges. 

Weak coupling (WC) g —> 0: Feynman perturbation theory is applicable in this case and so we use the Feynman 
diagram expansion to find 



^Pbuik = pk B T 1 - -g 
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The leading term is the free gas term and the second-order term is the familiar Debye-Huckel result in its ID variant. 
Note that there is no 0(g 2 ) term; this is cancelled by the counter term in Z(g). 

The strong and weak coupling dependencies of the bulk pressure Pbuik on g compare well with the exact solution 
of the problem. Both approximations are accurate across a wide range of g in their regime of validity. More details 
can be found in [4]. 



B. Finite L: exact methods 



For finite L we expand the kernel K((/)q, c/)l; L) over periodic eigenfunctions of the Mathieu equation. We can then 
use a numerical approach for eigenfunctions/eigenenergies which will give an exact solution for all L. This method is 
described fully in [4] and we do not delve into details here. It gives the same answers as the Fourier approach that 
we describe below. 

The Fourier method for obtaining an exact solution to problems in ID is more general than the Schrodinger approach 
since it works also when the Hamiltonian is not hermitian, which is the case for the counterion gas considered in the 
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next section. It also forms the basis for the transfer matrix method. The theory is periodic under <fi — >• <fi + 2ir and 
we can define 

/ n=oo 
K(4>,<f>';x)f(4>')d4>' = e x / 2 9 £ b n (x)e in <^, 
n= — oo 

where the coefficients b n (x) obey the evolution equation 

^ = -n 2 b n + -^On+i + 6 n -i - 26 n ). 

This is the Fourier version of the Schrodinger equation but can be derived generally from the convolution property of 
the Schrodinger kernel. The partition function can then be obtained from 

r 2ir 2 n=oo . n 

The exact solution for the disjoining pressure as a function of the separation L for different values of the surface 
potential strength parameter A clearly predicts a collapse transition to a Newton black film that can not be accounted 
for by the mean-field theory, which we address next. 



C. Classical or mean-field (MF) theory 



Standard variational methods applied to the expression for the partition function gives the classical MF equation: 
the Poisson-Boltzmann (PB) equation for (j) c i(x) as the saddle point equation of the corresponding field theory. In 
this case the disjoining pressure P is given by the value of ion density at the midpoint x = L/2 between the bounding 
surfaces. The MF theory predicts that universally P > 0, contrary to our exact result and also to experiment; it does 
not predict any collapse transition, which is thus obviously a consequence of the non-MF correlation effects and is 
intrinsically a fluctuation phenomenon. 



IV. COUNTERIONS BETWEEN CHARGED SURFACES 



The ID model here is a Coulomb gas of counterions confined between two oppositely charged surfaces; the system 
is overall neutral. We compare exact results with strong and weak coupling calculations, which are the same as in a 
3D system. More details can be found in |7]. 

The system is shown in Fig. [I] (middle) and consists of N counterions, each of valency q, with surface charges o\ 
and (72, respectively. We define £ = O2/01, with — 1 < £ < 1, and define a = 1/(1 + Q. The ID Bjerrum length is 
Ib = 2kBTe/e 2 1 and the Gouy-Chapman length is \i = n\ = Zse/g|(7i|, where we have chosen g\ to be non-zero and 
have H2 = /VICI- The electrostatic coupling constant, g, is then given by 
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where N — » oo corresponds to the MF/PB theory and N — >• 1 to the SC theory. The partition function is derived as 
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The integral over d(j)(0) ensures charge neutrality: Nqe + <ji +<J2 = 0. The corresponding Hamiltonian (Feynman-Kac) 
and the partition function are then 
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A. Exact results 



In this system H is not hermitian because the counterions are, by definition, of one charge only. We therefore 
analyze the model using the Fourier method. We exploit periodicity in H of <\> — >> <p + 2tt in order to write 

f(<l>; x) = e- xH f(4>; 0) with f(<f>; x) = £ b(n, x)e in ^ . 

By introducing g\ = — Miqe, a 2 = —M 2 qe^ Mi = Int(aiV), M2 = N — Mi — 1, r/i = aiV — Mi, and 772 = 1 — tji an d 
a = 1/(1 + C), we can derive the Fourier evolution equation from surface 2 to surface 1 in the form 

dHn ^ L) = -^^^e%(n,M 2 ,L) + K n-l,M 2 ,L), (4) 

with 6(n,M2,0) = 5 n ,-M 2 - This Fourier evolution equation can be integrated numerically and the corresponding 
partition function 

Z((7i,C,JV,L)=6(M 1 + l,M 2 ,L), 

and disjoining pressure 

P m (Mi + 1 - r\i) 2 2 2 b(M u M 2 ,L) 

can be evaluated exactly. Since the second term in the above equation can be seen to be just the counterion density 
at the boundary of the system, the above form of the pressure is thus a clear example of the contact value theorem; 
it connects the pressure with the value of the particle density at the confining wall of the system. 



B. Weak coupling 

We consider the WC expansion g — >> which is equivalent in the lowest order to the MF/PB theory. In the d > 1 
case, the MF theory treats the potential field </>(x) as constant in the directions transverse to the normal to the 
bounding interfaces, and so the results are independent of the dimensionality. 

The leading contribution arises from the saddle-point configuration c/)q(x) = ii^o(x) with real. The PB equation 
and the boundary conditions have the form 

dx z 



with 



dx 



dx 



cr 2 Pqe. 



o 

The leading PB contribution to the disjoining pressure, P , is then expressed as 

^-vW(f) 2 + -«- 

where po(x) is the density of counterions between the boundaries, given by the standard Boltzmann form po(x) = 

Ce~^^ x \ where C is a normalization constant. This furthermore implies that the MF/PB disjoining pressure P is 
obtained as follows: When the pressure is repulsive (P > 0), we have P = n 2 (j\Y 2 /2, where T satisfies 

tan(rL) = p^^, 

and when the pressure is attractive (P < 0), which may be the case within the MF/PB theory only for £ < 0, we 
have P = — ji 2 (j 2 Y 2 /2, where V is now given as a solution of 
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FIG. 2: Rescaled disjoining pressure, 2P/al, for a ID counterion gas between charg ed su rfaces as function of the rescaled 
intersurf ace se paration L/fi. Thick (red) solid lines represent the PB re sult (Section |IVB| ), dashed lines are the SC results 
(Section IV C) and thin solid lines are the exact results (Section IV A) compared with MC simulations data (symbols) at 
different numbers of counterions N and for £ = 0.5 (left pane) and £ = —0.5 (right pane) [7]. 



C. Strong coupling 

The strong coupling limit is formally identical to the one-particle limit pQ . In the present case it is easily evaluated 
from the partition function in the case of a single counterion in the system. The partition function in an explicit 
one-particle form leads to the disjoining pressure 

^) = ^(-^(l + C 2 ) + ^(l-C 2 )coth 

The range of validity of this limiting expression is of course defined by the number of counterions in the system. As 
this number decreases towards one, N —> 1, the above expression for the disjoining pressure becomes exact. 
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D. Comparison 



Both the weak and strong coupling approximations are independent of dimension d and the comparison with the 
exact results can test their validity. For symmetric surface charges (£ = 1) the PB/MF pressure is positive (repulsive) 
for all intersurface separations, whereas the SC expansion and the exact result for N = 1 predict attraction at large 
separations; this distinction holds for < £ < 1. For the asymmetric configuration with £ < 0, there is little difference 
between the different approaches; on trivial grounds there is attraction for large separations but there is repulsion for 
sufficiently small separations, see Fig. [2j where a comparison is made with Monte- Carlo (MC) simulations at different 
numbers of counterions N [7]. 



V. IONIC LIQUID LATTICE CAPACITOR 

In the models above the ions have been chosen to be point-like. Here we address the question of changes wrought 
by their finite size. In this case the system consists of a ID lattice of M sites with spacing a, with the i— th site, 
< i < M, occupied by ion with charge qSi with Si G [— g, 0,g], see Fig. [I] (bottom). Within this model the finite 
ion size is ~ a, which is crucial to the phenomena observed in experiments on confined ionic liquids. 

The configuration described is one of the ID ionic liquid capacitor. The external fields are imposed either by fixing 
the charges of the boundaries at i = — 1 and i = M to be respectively, or by imposing a fixed voltage/potential 

difference, Av, across the capacitor. More details can be found in [8]. 

The electrostatic Hamiltonian in this case is expressed through a spin-like variable Si = 0, ±1 

m = -ff\i-j\s iSj 7 = ^. 

i.j=0 
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After a Hubbard- Stratonovich transformation this yields the action 

M-2 , , _ , , 2 M ~ 1 

j=o 7 i=o 

The system includes boundary charges ±gQ a ^ sites —1, M. The electrostatic potential is defined as V = —i(j)/f3q. In 
limit a — > 0, g/a fixed, the MF equations obtained from the saddle-point of the above field action reduce to those of 
Kornyshev [9 and Borukhov et al. [TO] . 

For non-zero a the action is not positive definite for \i > 0.5 and so we seem to have a sign problem and certainly 
cannot use the Schrodinger approach a priori. Nevertheless, in the case of ID the partition function can be computed 
exactly by using the transfer matrix approach, with the Fourier method described earlier. This can be seen as follows: 
write yi = fa and define 

p 1/s W) = ^=e-^-y'?i\ 

V^7 



with 

K(y,y') = J dzp l l 2 (y,z)[l + n™s(z)}p l l 2 {z,y') (5) 
with Kf(y) = dy' K(y J y , )f(y / ). The free energy for the fixed Q ensemble, ^q, then follows as 

-pa Q= r dx e iQ* Li/2^1/21 e -i Qx = {i/, q \K m \i/>q) 

J —TV 



e 



with 

(iP Q \y) =e«lv, (y\K\y') = K(y,y') 
The conjugate free energy for the fixed Av ensemble, Qav, then follows from a Legendre transform, 



while the capacitance Cav is obtained from the first derivative of 8{Q)av w.r.t. Av. and can thus be calculated 
directly from the partition function. 



A. Results 



The transfer matrix and Fourier approach can be formulated in order to evaluate the free energy explicitly. Details 
of this procedure can be found in Ref. [8 . Enthalpy Qm — &m + ^ftuik, the disjoining pressure P = Qm — Gm+i, 
and the capacitance Cav, can all be calculated as a function of /i, Q, Av. 

We show explicitly only the capacitance results, Cav, as a function of Av in Fig. [3] both for large ji and small 
li. For large fi the curve shows shows the typical "bell" shape in contrast to the curve for smaller /i, which shows 
the non-monotonic "camel" shape and so Cav nas a minimum at the point of zero charge confirming the Fermi MF 
results of Kornyshev [9] . 

For smaller 7 (increasing T) the periodic non-monotonicity both for large /i and small fi disappears and the solution 
approaches the Fermi MF result of Kornyshev [9 . It is interesting that the exact solution dances around the Kornyshev 
solution with an ever increasing amplitude but the system nevertheless always remains thermodynamically stable, as 
can be straightforwardly ascertained. 



VI. LESSONS 



We have demonstrated that in ID one can use the Schrodinger approach for continuum models of Coulomb fluids, but 
that for discrete models a more general approach is needed which exploits the transfer matrix and the periodicity of the 
field to use Fourier methods. We tested the PB/MF and the strong coupling limiting expressions and demonstrated 
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FIG. 3: (Left) Bell shaped Ca v as a function of Av for large \i — 0.5. (Right) Camel shaped Ca v as a function of Av for small 
[i — 0.03. In both cases, we have 7=1. The dashed line is the MF result [8]. 



that they need correcting although the exact analytic result clearly supports the two limiting analyses. We also 
confirmed that the MF theory does not capture the important effects which are due to correlations, either the 
attractive intersurface forces in the case of a counterion-only system or non-monotonic periodic variation of the 
capacitance in the confined ionic liquid case. 
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